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In this article we consider a new derivation mode for generalized communicating P systems (GCPS) 
corresponding to the functioning of population protocols (PP) and based on the sequential derivation 
mode and a fairness condition. We show that PP can be seen as a particular variant of GCPS. We 
also consider a particular stochastic evolution satisfying the fairness condition and obtain that it 
corresponds to the run of a Gillespie's SSA. This permits to further describe the dynamics of GCPS 
by a system of ODEs when the population size goes to the infinity. 

1 Introduction 

The notion of a generalized communicating P system was introduced in ll22l . with the aim of providing 
a common generalization of various purely communicating models in the framework of P systems. 

A generalized communicating P system, or a GCPS for short, corresponds to a hypergraph where 
each node is represented by a cell and each edge is represented by a rule. Every cell contains a multiset 
of objects which - by communication rules - may move between the cells. The form of a communi- 
cation rule is (a,i)(b,j) — > (a,k)(b,l) where a and b are objects and i,j,k,l are labels identifying the 
input and the output cells. Such a rule means that an object a from cell i and an object b from cell j 
move synchronously to cell k and cell I, respectively. In this respect, the model resembles the Petri Net 
formalism [18] where tokens from various input places come along together to fire a given transition and 
then fork out to destination places, see ll22l |4 ] for more details. 

Depending on the communication rules form, several restrictions on communication rules (modulo 
symmetry) can be introduced. Due to the simplicity of their rules, the generative power of such restricted 
systems is of particular interest and it has been studied in detail. In ll22l [U |71 it was proved that 
eight of the possible nine restricted variants (with respect to the form of rules) are able to generate any 
recursively enumerable set of numbers; in the ninth case only finite sets of singletons can be obtained. 
Furthermore, these systems even with relatively small numbers of cells and simple underlying (hyper- 
graph) architectures are able to achieve this generative power. In Q a further restriction is introduced 
by considering that the alphabet of objects is a singleton (like in Petri Nets) and it is shown that the 
computational completeness can be achieved in four of the restricted variants. 

Population protocols (PP) have been introduced in flU (see for a survey) as a model of sensor net- 
works consisting of very limited mobile agents with no control over their own movement. A population 
protocol corresponds to a collection of anonymous agents, modelled by finite automata, that interact with 
one another to carry out computations, by updating their states, using some rules. Their computational 
power has been investigated under several hypotheses in most of the cases restricted to finite size popula- 
tions. In particular, predicates stably computable in the original model have been characterized as those 
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definable in Presburger arithmetic. The article 0] studies the convergence of PP when the population 
size goes to the infinity. 

The evolution of a PP follows a particular fairness condition: an execution is fair if for all configu- 
rations C that appear infinitely often in the execution, if C is predecessor of a configuration C , then C 
appears infinitely often in the execution. We consider such a condition in the case of GCPS systems 
and obtain a new derivation mode which we call fair sequential mode (fs-mode). We further study the 
dynamic behaviour of the system in this mode. Among several possible evolution strategies we con- 
sider a stochastic strategy that satisfies the fairness condition and we obtain that the evolution of the 
system corresponds to a run of the Gillespie stochastic simulation algorithm (SSA). Using the corre- 
spondence between SSA and ODEs (assuming mass-action kinetics) we show that the dynamics of the 
system can be represented by a system of ODEs when the population size goes to the infinity. We also 
consider the converse problem and we give sufficient conditions for a system of ODEs to be represented 
by a GCPS system working in concentration-depended stochastic implementation of the fs-mode. We 
consider several examples of GCPS simulating Lotka-Volterra (predator-prey) behaviour or computing 
approximations of algebraic numbers. 

2 Background 

In this section we recall some basic notions and notations used in membrane computing, formal language 
theory and computability theory. For further details and information the reader is referred to |[T6l[T7l[T9l . 

An alphabet is a finite non-empty set of symbols. For an alphabet V, we denote by V* the set of all 
strings over V, including the empty string, X. The length of the string x G V* is the number of symbols 
which appear in x and it is denoted by \x\. The number of occurrences of a symbol a G V in x G V* is 
denoted by \x\ a . If x G V* and U C V, then we denote by \x\u the number of occurrences of symbols 
from U in x. 

A finite multiset over V is a mapping M : V — > N; M(a) is said to be the multiplicity of a in M (N 
denotes the set of non-negative integers.) A finite multiset M over an alphabet V can be represented by 
all permutations of a string x = af^^a^"^ . . . a£ G V*, where aj G V, 1 < j < n; x represents M 
in V*. If no confusion arises, we also may use the customary set notation for denoting multisets. The 
size of a finite multiset M, represented by x G V* is defined as Z ae y |x| a . 

2.1 P Systems 

Next we recall the basic definitions concerning generalized communicating P systems l22l . 

Definition 1. A generalized communicating P system (a GCPS) of degree n, where n > 1, is an (n + 4)- 
tuple n = (0,E,w\, . . .,w n ,R,h) where 

1. O is an alphabet, called the set of objects of IT; 

2. £CO; called the set of environmental objects of IT; 

3. Wi G O*, 1 < i < n, is the multiset of objects initially associated with cell i; 

4. R is a finite set of interaction rules (or communication rules) of the form (a,i)(b,j) — > (a,k)(b,l), 
where a, b G O, < i, j, k,l <n, and if i = and j = 0, then {a, b} C\(0\E) 7^ 0; i.e., a ^ E and/or 
b^E; 

5. h G {1, ... ,n} is the output cell. 
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The system consists of n cells, labelled by natural numbers from 1 to n, which contain multisets of 
objects over O; initially cell i contains multiset w, (the initial contents of cell i is w,-). We distinguish an 
additional special cell, labelled by 0, called the environment. The environment contains objects of E in 
an infinite number of copies. 

The cells interact by means of the rules (a,i)(b,j) — > (a,k)(b,l), with a,b £ O and < i,j,k,l < n. 
As the result of the application of the rule, object a moves from cell i to cell k and b moves from cell j 
to cell /. If two objects from the environment move to some other cell or cells, then at least one of them 
must not appear in the environment in an infinite number of copies. Otherwise, an infinite number of 
objects can be imported in the system in one step. 

A configuration of a GCPS IT, as above, is an (n + l)-tuple (zo,zi, . . . ,z n ) with zo £ (0\E)* and zi £ 
O* , for all 1 < i < n; zq is the multiset of objects present in the environment in a finite number of copies, 
whereas, for all 1 < i < n, Zi is the multiset of objects present inside cell i. The initial configuration ofYl 
is the tuple (X,Wi, . . ■ ,w„). 

Given a multiset of rules ffl over R and a configuration u = (zo,Zi,- ■ ■ ,z n ) of IT, we say that & is 
applicable to u if all its elements can be applied simultaneously to the objects of multisets zo,Zi, ■■■,z n 
such that every object is used by at most one rule. Then, for a configuration u = (zo,Zi, ■ ■ ■ ,z n ) of n, 
a new configuration u' = (Zq,z\, ■ ■ ■ ,z'„) is obtained by applying the rules of R in a non-deterministic 
maximally parallel manner: taking an applicable multiset of rules M over R such that the application 
of ffl results in configuration u' = (z' ,z\, ■ ■ ■ ,z' n ) and there is no other applicable multiset of rules S$ 
over R which properly contains M. 

It is also possible to replace the maximally parallel strategy of rule application by other strategies, 
called derivation modes (in the context of the present paper, the terms mode and strategy are used in- 
differently). A derivation mode lies in the heart of the semantics of P systems and it permits to specify 
which multiset among different possible applicable multisets of rules can be applied. When P systems 
were introduced, only the maximally parallel derivation mode was considered which states that corre- 
sponding multisets should be maximal, i.e., non-extensible. With the appearance of the minimal parallel 
derivation mode [6] the concept of the derivation mode had to be precisely defined and [ 10] presents a 
framework that permits to easily define different derivation modes. 

One application of a multiset of rules satisfying the conditions of a derivation mode represents a 
transition in IT from configuration u to configuration u'. A transition sequence is said to be a successful 
generation by IT if it starts with the initial configuration of IT and ends with a halting configuration, i.e., 
with a configuration where no further transition step can be performed. 

We say that IT generates a non-negative integer n if there is a successful generation by IT such that n 
is the size of the multiset of objects present inside the output cell in the halting configuration. The set of 
non-negative integers generated by a GCPS IT in this way is denoted by N(IY). It is also possible to use 
GCPS as acceptors, in this case an input multiset is accepted if the system halts on it. 

In [22] it is shown that GCPS are able to generate all recursively enumerable languages. Moreover 
this result can be obtained by using various restrictions on the type of rules (i.e. induced hypergraph 
structures), on the number of membranes and on the cardinality of the alphabet. We refer to 11221 1811711 for 
more details. 

If the cardinality of the alphabet O is equal to one, then we refer to the corresponding symbol as a 
token (denoted by •). Hence, we assume that O = {•}. We observe that such systems are similar to Petri 
Nets having a restricted topology. This is especially visible if a graphical notation is used. However, the 
maximal parallelism and the concept of the environment are specific to P systems, so we place this study 
in the latter framework. A converse study of P systems from the point of view of Petri Nets can be found 
in ifTTi . For more details on Petri Nets and membrane computing we also refer to ifTTl . 



85 



In this article we shall consider the dynamics of the configuration of GCPS, so we are no more 
interested in computation (and halting evolutions). 

2.2 Population Protocols 

We give below the definition as it appears in 10. A protocol is given by (<2,r,i, ft), 5) with the following 
components. Q is a finite set of states. £ is a finite set of input symbols. 1 : £ — > Q is the initial state 
mapping, and ft) : Q — > {0, 1} is the individual output function. 8 C Q 4 is a joint transition relation that 
describes how pairs of agents can interact. Relation 8 is sometimes described by listing all possible 
interactions using the notation (^1,^2) — > (?'i>?2)' or even tne notation q\q2 — > q[<l2> f° r (?i > I2i 4\ > 4%) £ 
8 (with the convention that (91,^2) — > {qiiQi) when no rule is specified with (^1,172) hi the left hand 
side.) 

Computations of a protocol proceed in the following way. The computation takes place among n 
agents, where n > 2. A configuration of the system can be described by a vector of all the agent's states. 
The state of each agent is an element of Q. Because agents with the same states are indistinguishable, 
each configuration can be summarized as an unordered multiset of states, and hence of elements of Q. 

Each agent is given initially some input value from E: each agent's initial state is determined by 
applying 1 to its input value. This determines the initial configuration of the population. 

An execution of a protocol proceeds from the initial configuration by interactions between pairs of 
agents. Suppose that two agents in state qi and q2 meet and have an interaction. They can change into 
state q\ and q' 2 if {q\ , q2, q\ , q' 2 ) is in the transition relation 8. If C and C are two configurations, we write 
C — > C' if C' can be obtained from C by a single interaction of two agents: this means that C contains two 
states q\ and q2 and C is obtained by replacing q\ and q2 by q\ and q' 2 in C, where (^1,^2,^11^2) e ^- An 
execution of the protocol is a (potentially infinite) sequence of configurations Co,C\,C2, ■ ■ ., where Co is 
an initial configuration and C, — > C,+i for all i > 0. An execution is fair if for all configurations C that 
appears infinitely often in the execution, if C — > C for some configuration C', then C' appears infinitely 
often in the execution. 

At any point during an execution, each agent's state determines its output at that time. If the agent 
is in state q, its output value is 0){q). The configuration output is (respectively 1) if all the individual 
outputs are (respectively 1). If the individual outputs are mixed 0s and Is then the output of the 
configuration is undefined. 

Let p be a predicate over multisets of elements of Z. Predicate p can be considered as a function 
whose range is {0, 1} and whose domain is the collection of these multisets. The predicate is said to 
be computed by the protocol if, for every multiset /, and every fair execution that starts from the initial 
configuration corresponding to /, the output value of every agent eventually stabilizes to p(I). 

The following was proved in HUH: 

Theorem 1 (UJ |3). A predicate is computable in the population protocol model if and only if it is 
semilinear. 

Recall that semilinear sets are known to correspond to predicates on counts of input agents definable 
in first-order Presburger arithmetic IPT31 . 

2.3 Gillespie Algorithm 

A usual abstraction in the simulation of biochemical systems consists in considering the system (e.g., a 
bacterium) as a homogeneous chemical solution where the reactions of the model are taking place. D.T 
Gillespie has proposed in lfl2l an algorithm for producing the trajectories of such a chemical system 
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by computing the next reaction and the elapsed time since last reaction occurred. Let fx be a chemical 
reaction. The probability that /I takes place during an infinitesimal time step is proportional to: 

• Cp, the stochastic reaction constan^of reaction /i; 

• h s ^ , the number of distinct molecular combinations that can activate reaction /I ; it depends on the 
current chemical state 5; 

• dz, the length of the time interval. 

Gillespie proved that the probability P(z,fl\S)dz that, being in a chemical state 5, the next reaction will 
be of type pt and will occur in the time interval (t + z, t + z + dz) is: 

P(z,ix\S)dz = af l e- a o' t dz 

where a s ^ = c^h^ is called the propensity of reaction ji., and a$ = Y*v a v i s tne combined propensity of 
all reactions. 

This probability leads to the first straightforward Gillespie's exact stochastic simulation algorithm 
(SSA) called theirs? reaction method. From a chemical state S, it consists in choosing an elapsed time 
z for each reaction /I according to the probability P (z,/J.\S). The reaction with the lowest elapsed time 
is selected and applied on the system making its state evolve. A new probability distribution is then 
computed for this new state and the process is iterated. 

The Gillespie's SSA gives a way to simulate a continuous-time Markov chain with the states cor- 
responding to the states of the system and with transitions between states corresponding to a single 
occurrence of a reaction. The probability for a transition between two states 5 and S' corresponding to 
the application of rule pt is defined as a^/a^. In the following, we drop the mention of the current state 
S in these notations. 

3 Fair Sequential Derivation Mode 

In this section we are interested in the relation between PP and GCPS. We show that in terms of structure 
PP and GCPS are quite similar, the main differences concern the environment and the derivation mode. 
We define a new fair sequential mode for GCPS and hence we are able to encode any PP in a GCPS w.r.t. 
their dynamics. We then remark that GCPS with stochastic and Gillespie-like strategies are part of this 
new class of GCPS and we propose their use for simulations of population behaviours. 

It can be easily seen that both PP and GCPS are particular instances of multiset rewriting. Indeed, 
in both cases the underlying data structure is multiset (obtained in a direct way for PP and by attaching 
the indices of membranes to the objects in GCPS) and the evolution rules clearly correspond to multiset 
rewriting rules with both left hand and right hand sides of size two. So, the translation of a PP to a 
one-symbol GCPS can be easily done as follows. Given a PP with set of states Q (for convenience we 
suppose that Q = {1, . . . ,«}) and transition relation 8 in an initial configuration Co, the corresponding 
GCPS n = (0,E,w\, . . .,w„,R, 1) can be defined as: 

• = £ = {•}, 

• w q = » k ,k = \C \ q for any q G Q, 

• r = {(•,<?!)(•, <? 2 ) -> {•,q\){»,q' 2 ) I qm -> qWi e 

'Evaluating the stochastic constants is one of the key issues in stochastic simulations of biochemical reactions. 
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The above system encodes each state q of PP by a token • present in membrane labelled by q. Since we 
are interested in the dynamics of the system, no output membrane is necessary. The above construction 
covers the core of PP and to obtain the complete equivalence encoding and decoding functions 1 and CO 
shall be used in the same way. 

It is not possible to do a similar encoding of GCPS with PP because PP always deal with finite 
multisets and GCPS can use the infinite multiset corresponding to the environment. However, any GCPS 
having no rule involving the environment can be translated to PP in a similar way. 

We remark that the biggest difference between PP and GCPS is given by the evolution step, i.e., 
by the derivation mode. For GCPS, mainly the maximally parallel derivation mode is investigated with 
several attempts to investigate asynchronous or minimally parallel derivation mode, see [ 17"] for more 
details. The derivation mode of PP is very particular - it corresponds to a sequential strategy where only 
one rule is applied at each step, like in Petri Nets, but with an additional fairness condition. 

We can consider such a strategy in GCPS case as well. More precisely we consider fair computations: 
a GCPS computation is fair, if for any configuration u that appears infinitely often in the computation, 
then any configuration u' , such that u => u' in sequential application, also appears infinitely often. We 
shall call such computational strategy a. fair sequential derivation mode (shortly fs-mode). 

From these considerations, it is trivial to observe that PP are similar to GCPS in fs-mode with only 
one symbol in the alphabet. If we consider an encoding function 1 like for PP and the halting condition 
corresponding to the stabilization of the ft)-image of the configuration, then as an immediate consequence 
of HI HI, we obtain that any GCPS working in fs-mode and that does not have any rule involving the 
environment can only accept semilinear sets. 

Conversely, we also obtain that any PP working in maximally parallel mode {i.e., a maximally parallel 
number of interactions can happen at each step) are computationally complete if the number agents in 
some particular state q$ is going to the infinity. 

From now on, we only speak of PP in terms of their associated one-symbol GCPS in fs-mode. 

3.1 FS-Mode and Stochastic Evolution 

Although powerful the definition of the fairness remains obscure. Let try to clarify it. When assuming 
that the number of configurations is finite (this is the case for classical PP for example), the definition 
can be easily rephrased as follows: a computation uq => u\ =>■ . . . is fair if 

• there exists a non-negative integer N such that configuration belongs to a terminal strongly 
connected component of the state graprQ and 

• any state of this terminal strongly connected component appears infinitely often in the execution. 

There are many possible evolutions of the system in the fs-mode. One example of such an evolution is 
to choose at each step a rule that leads to a configuration that either never was visited previously or was 
not visited for some time greater than k, k > (if possible). 

Among all possible evolutions, Markovian processes feature prominently since they respect the fair- 
ness condition (it is well-known that Markovian processes leave non-terminal strongly connected com- 
ponents with probability 1), they do not require any history or global knowledge on the state space, and 
they provide a modelling tool useful in many domains (like in the simulation of population behaviours 
or in distributed algorithmics). Such a Markovian process corresponds to a labeling of each state graph 
arrow u ^> p u' by a static probability p that only depends on configuration u. Here are two examples of 
such Markovian processes: 

2 In this directed graph, nodes correspond to the configurations u, and two nodes u and u' are directly linked if u =>■ u'. 
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1. Equiprobable evolutions: u =^V*» u' where k u denotes the cardinality of the set {u' \ u =>• u'}. 

2. Concentration-dependent evolutions: u =$- Pr u' where r denotes the applied rule and p r is propor- 
tional to h r , the number of distinct combinations of tokens that activate r, with a proportionality 
coefficient that only depends on r. Assuming that r = (*,qi)(*,q 2 ) — > (*,q'i){*i a 2)> the number 
h r is given by 

|"M"|? 2 if <?i ^°!?2 /0,<7i ^q 2 

\ U \qA\ U \qi ~ 1) i f <?1 ^ 0,^1 = <?2 



ft, 



191 

1 92 



if q 2 
if 91 








(1) 



The two last cases hold when the environment (containing an infinite number of tokens) is involved 
in the rule. 



3.2 FS-Mode GCPS modelling Population Dynamics 

Assuming that for a given rule proportionality coefficients are the same for all configurations, the con- 
centration-dependent strategy directly corresponds to a run of the Gillespie's SSA. Thus, we advocate 
that GCPS in fs-mode provide a good theoretical tool for studying population behaviours. 

A paradigmatic example illustrating how GCPS allows a well suited specification of population be- 
haviours consists of the description of a process inspired by the Lotka-Volterra model. 



The Lotka-Volterra Model. The Lotka-Volterra process was introduced by Lotka as a model of cou- 
pled auto-catalytic chemical reactions, and was investigated by Volterra as a model for studying an 
ecosystem of predators and preys (9). This model specifies how two coupled populations (of chemi- 
cals or individuals) Y\ (the preys) and Y 2 (the predators) behave. In lfT2l . D.T Gillespie proposes the 
study of this system derived from the following ODEs 

d ^ = {cx-c 2 Y 2 )Y 1 ^ =(c 2 Yi-c 3 )Y 2 (2) 

Equivalently, the following chemical reactions 

Cl c 2 C3 

Y\ >2Y 1 Y\ +Y 2 >2Y 2 Y 2 * . (3) 



specify a model whose behaviour is described by ODEs system ([2]). The dynamics of these reactions 
is conveniently characterized using the predator-prey interpretation. The first rule states that a prey Y\ 
reproduces. The second rule states that a predator Y 2 reproduces after feeding on prey Y\ . Finally, the last 
rule specifies that predators Y 2 die of natural causes. Coefficients c, are the rates of the three reactions. 
The correspondence between the two models relies in the fact that the trajectories of the Gillespie's SSA 
tend to the solutions of the ODEs system given by the law of mass action on the reactions. This result is 
due to the particular application of the Kurtz's theorem iTPl to chemical systems. 



Lotka-Volterra GCPS Definition. The model above does not fill GCPS requirements since the first 
and last reactions are not pairwise interactions. We propose to extend reactions ([3]) by considering a 
renewable resource X for Y\ as a third specie^] the molecular level of X remains constant whatever is 

3 We use the same notation as in 1121 to express that the food resource X is assumed renewable. 
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Figure 1: Runs of the Lotka-Volterra model with renewable (left column) and not-renewable (right column) 
resources for different initial states (kinetics rates equal 1): Y\ = 500, Yi = 1500, X = X = 1000 (first row), Y\ = 
Y 2 = X = X = 1000 (second row), Y\ = Y 2 = 1000, J = X = 10000 (third row), Y\ = 1000, Y 2 = 0, X = X = 10000 
(fourth row). The solid red line represents preys, the dashed green line predators, and the blue dotted line resources. 
The two first rows show that both dynamics exhibit the same properties as presented in lfl2l (particularly, in second 
row, oscillations raise from an equilibrium initial state for the ODEs). The third row shows the difference in the 
dynamics when the resource size is ten times larger tharPGie population size. The last row shows the difference in 
the dynamics when the predator population is empty. The simulations have been done using the general simulation 
language MGS (http://mgs.spatial-computing.org) that allows an easy implementation of all models of 
the present article 1201 IT51 12T1 . 



its production or its consumption. The extended system of reactions is: 

C y ^2 C3 

X + Yi >2Y l Y1+Y2 >2Y 2 Y 2 +X >2X ( 4 ) 

The use of a pairwise interaction in the last reaction can be interpreted as a competition between the two 
predator behaviours: a predator in presence of preys eates and reproduces (second reaction); a predator 
in absence of prey (represented by the grass) dies (third reaction). Moreover, with the hypothesis that 
the number of X remains constant, the behaviour of this system is exactly described by ODEs ^ with 
ci = Xc\, c 2 = c' 2 and C3 = Xc' 3 . Thus, systems ^ and Q are equivalent in terms of dynamics. 

System ([4]) is only composed of pairwise interactions that satisfy condition 4 of Definition [T] Thus, 
it can be easily translated to a one-symbol GCPS, denoted Ulv, working in fs-mode (in concentration- 
dependent implementation) with rules R: 

(.,0)(.,1)-K.,1)(.,1) («,1)(«,0) (•,!)(•,!) 

(.,1)(.,2) -> (.,2)(.,2) (-,2)(.,1) -> (.,2)(.,2) 

(.,2)(.,0) -> (.,0)(.,0) (-,0)(.,2) -> (.,0)(.,0) 

where membrane indices 0, 1 and 2 represent the environment (an infinite source of X), the preys Y\ and 
the predators Y 2 , respectively. 

Let now consider the previously defined concentration-dependent evolution with probabilities = 
a^/ao for each pt G R with the propensity function =c^h^: is the rate of the corresponding reaction 
in ([3]> and is given by equation ([T]) accordingly to Ulv- The reader is invited to pay attention that even 
if the environment is an infinite source of X (instead of a constant one), the dynamics are well taken into 
account: rules involving the environment have probabilities that do not depend on the environment size, 
see equation ([!]). For example, the propensity of the first reaction is given by a\ =c\h\= c\ Y\ = c\XY\ 
as expected w.r.t. reactions ([4]). In this respect, any computation of n^y represents a run of the Gillespie's 
SSA of reactions Q. As a consequence, Ulv is an exact model of the original Lotka-Volterra system. 

It has to be remarked that Ulv cannot be described by any PP since the environment objects are 
involved in its definition. A possible specification of the Lotka-Volterra equations may be obtained 
within a PP by considering X as a not-renewable resource. Such a definition has been realized (taking 
reactions Q and substituting X by X.) However, due to the limitation of resource, this system does not 
respect the dynamics of equation ([2]) anymore. For example, without any predators, a population of preys 
stabilizes in this model, while in the original model it grows exponentially. Figure[T]gives some examples 
of simulations of the Lotka-Volterra model considering renewable and not-renewable resources. 



General Population Dynamics. It is possible to reverse the above method and to give a GCPS system 
whose population dynamics will correspond to some dynamics given by a system of ODEs, under the 
following conditions. Let us consider the ODEs system defined on set of variables {Y\, . . -,Yn} of the 
form 

dY 

-j7 =L«W-I(^7 + ^m (5) 

m j,k j 

where coefficients a l - k and by satisfy the following conditions: 

1. for all i,j,k, a'j k > and by > 0; 

2. for all j,k such that bjt 7^ 0, there exists either one index io such that a\ = 2b jt, or two distinct 
indices i\ and i 2 such that a\ = aj k = bj^; for any other index i, i ^ io or i ^ i\ and i 7^ i 2 , a l - k = 0. 
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The above conditions are sufficient to ensure that = 0. Then there exists a concentration-dependent 
fs-mode GCPS without rules involving the environment (i.e., a PP) whose behaviour is exactly described 
by ODEs Q when the population size goes to the infinity. Indeed, these equations correspond to the 
mass-action law of a set of rules such that for any j,k with bjk 7^ 

(•J)(*,k) (•,k)(*,k) or (;j)(;k) -+ b * (•,h)(*,h) 

according to the considered possibility of the above condition 2. The reader is invited to pay attention that 
these equations correspond to a wider range of dynamics than the dynamics of second-order chemical 
reactions with two products since they allow the specification of ordered interactions (e.g., involving a 
sender and a receiver as considered in the PP literature). This property also holds in PP and suggests that 
equations (|5]> exactly describe PP dynamics when the size of the populations tends to the infinity. 

It is obvious that a more general class of population behaviours is captured by concentration-dependent 
fs-mode GCPS since they have not to be conservative thanks to the environment. Following the idea of 
equivalence between systems ([3]) and ([4]) in terms of dynamics, equations Q can be extended with the in- 
troduction of a renewable variable Yq. This wider class of ODEs is supported by concentration-dependent 
fs-mode GCPS model. 

4 Computational Properties 

In this section, we focus on the original use of PP as a computational model of algebraic numbers 
proposed in (3). This article investigates the case where the computation is independent of the initial 
contents of the system. 

Using the GCPS terminology, the main idea of is to consider the result of a computation as a 
ratio between the number of tokens in certain membrane and the total number of tokens (without taking 
care of the environment) when the population size goes to the infinity and when the state of the system 
converges. The proposed work relies on the definition of a particular strategy of execution of the PP: 
a step of execution consists in sampling uniformly and independently of the past two distinct tokens in 
the membrane and let them interact in a sequential mode. This strategy is fair since it corresponds to a 
Markov process. The authors of the aforementioned article studied the Markov chain associated with PP 
and proved its equivalence to some system of ODEs at the limit. 

We remark that the same kind of result directly arises from considerations of Section |3.2| since this 
computational model is captured by one-symbol GCPS working in fs-mode with Gillespie concentration- 
dependent implementation. Indeed, the above execution strategy exactly corresponds to a Gillespie's 
SSA run where the stochastic constants equal 1 for all rules. Thus, the study of the model corresponds 
to the investigation of the sensibility of the associated ODEs system. Let us illustrate this point by 
considering the running example of 

(;p){;p) -> (•,p)(;m) 

(;p)(;m) -> (;p)(;p) 
{;m)(;p) ->• (;p)(;p) 
(•,m)(»,m) -> (;p)(;m) 

where symbols p and m identify two membranes. It has been shown that the ratio , where p (resp. m) 
is the size of the membrane p (resp. m), converges to -4= when the population size goes to the infinity. 
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Accordingly to equations ([5]>, we associate ODEs with this GCPS as follows 

dY ^-Y 2 + 2YY -Y 2 — --Y 2 -2YY+Y 2 

The stable states of this system are obtained when the two equations vanish, that is, when either Y m = 
— (y/2 + VjYp or Y m = (\/2 — l)Y p . The first solution is incoherent since it involves a negative size of 

Y 1 

population. The second solution trivially leads to the expected result Y £ Y = 

5 Conclusions 

In this article we investigated connections between population protocols and generalized communicating 
P systems. The two models share the same multiset structure and the same type of rules. Traditionally 
PP are used to study population dynamics in the context of distributed algorithmics while GCPS are 
investigated for the computational properties. 

By incorporating the derivation mode from PP into GCPS framework we obtained a strict inclusion of 
PP in GCPS working in fs-mode. We then took a particular implementation of the fs-mode corresponding 
to a run of the Gillespie's SSA and we obtained that the dynamics of the systems can be described by 
the corresponding system of differential equations. Different questions then could be explored, like the 
investigation of the conditions ensuring that the system reaches a stable state regardless of its initial state 
or ensuring that a stable state is never reached for any initial configuration. GCPS are in this sense easier 
to handle than PP because of the environment that permits to easily simulate the equivalent of creation 



or degradation reactions. Section 3.2 also considers the converse problem of the construction of a GCPS 
system exhibiting a particular behaviour given by a systems of ODEs. It would be interesting to see if 
the given sufficient conditions are also necessary. A mathematical challenge resulting from Section [4] 
is whether for any algebraic number x £ [0..1] there is a GCPS working in concentration-dependent 
evolution implementation of the fs-mode that converges to x. 

We remark that the presented results hold only in the concentration-dependent implementation of the 
fs-mode. By taking an equiprobable implementation the results are completely different. 

Since Petri Nets can be seen as multiset rewriting, it is clear that the results of this paper can be 
translated to this domain (for Petri Nets with specific type of rules and an additional fairness strategy). 

We think that the fs-mode has interesting properties that should be further explored. As showed in the 
article, the fairness condition is in some sense similar to a stochastic evolution, so it could be preferable 
to consider this condition instead of a stochastic behaviour. Another interesting property of the proposed 
stochastic implementation is that Gillespie's SSA introduces an explicit continuous time and discrete 
events in the model, which do not appear in a GCPS description. 
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